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We measure the sum of the neutrino particle masses using the three-dimensional galaxy 
power spectrum of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) Data Re- 
lease 9 (DR9) CMASS galaxy sample. Combined with the cosmic microwave background 
(CMB), supernova (SN) and additional baryonic acoustic oscillation (BAO) data, we find 
upper 95 percent confidence limits of the neutrino mass Hm^^ < 0.340 eV within a flat 
ACDM background, and Sm^y < 0.821 eV, assuming a more general background cosmo- 
\ logical model. The number of neutrino species is measured to be A^eff = 4.308 ± 0.794 

and A^eff = 4.032l^i^^ for these two cases respectively. We study and quantify the effect 
of several factors on the neutrino measurements, including the galaxy power spectrum bias 
model, the effect of redshift- space distortion, the cutoff scale of the power spectrum, and the 
choice of additional data. The impact of neutrinos with unknown masses on other cosmolog- 
ical parameter measurements is investigated. The fractional matter density and the Hubble 
parameter are measured to be = 0.2796 ± 0.0097, Hq = 69.72t^|? km/s/Mpc (flat 
ACDM) and flu = 0.279StoS[ll, Hq = 73.78t^;j^ km/s/Mpc (more general background 
model). Based on a Chevallier-Polarski-Linder (CPL) parametrisation of the equation-of- state 
w of dark energy, we find that w = —lis consistent with observations, even allowing for neu- 
trinos. Similarly, the curvature and the running of the spectral index as are both consistent 
with zero. The tensor- to- scalar ratio is constrained down to r < 0.198 (95 percent CL, flat 
ACDM) and r < 0.440 (95 percent CL, more general background model). 
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1 INTRODUCTION 

The detemiination of the neutrino mass is one of the most im- 
portant tasks in modem science. High-energy experiments probe 
the mass differences between neutrino species through neutrino 
oscillations. Latest measurements give squared mass differences 
for the solar neutrino and the atmospheric neutrinos of Am?2 = 



7.59: 



X 10"^ eV 



^ a nd Amjc. = 2.35 



+0.12 
n nc 



X 10"''eV\ re- 



spectively jSchwetz et al.ll201ll : lFodi et al.ll201ll 

While there is no way to measure the absolute neutrino mass 
in these oscillation experiments this can be measured using cosmo- 
logical observations. The redshift of the matter-radiation equality 
depends on the summed mass of neutrino particles, leaving an im- 
print on the cosmic microwave background (CMB) angular power 
spectrum, and the matter transfer function. Massive neutrinos also 
affect the cosmological growth factor, inhibiting growth as their 
velocity dispersion slows down their mean flow into struct ures. Re- 
views of th e effects of neutrin os in cosmology are given bv lPolgovl 
(l2002h and lLesgourgues & P astor (2006). 

Given the number of recent galaxy surveys, and different 
options for modelling non-linear effects in the measured galaxy 
power spectrum, it is not surprising that there have been a signifi- 
cant number of previous efforts to measure the neutrino mass using 
cosmological probes, {e.g., Elgar0v et al. 2002; Lewis & Bridle] 
2OO2I: lAllen et all l2003l: iHannestadl l2003l : ISoergel et all |200; 
Bargeretal.1 l2004l: ICrottv et all |2004 iHannestad & Raffelt 



2004 



2005 



2007 



Tegmark et aP 120041 : E'lgar0v & Lahavl l2005l: IHannestad 
iSeliak et al. 2005; Goobar et al. "200^: Spergel et al 
Jxiaetal. 2007; Komatsu et al. 2009; Tereno et al. 200< 
ISaitoetalJ I2OIII: ISwanson et all 



Reidetal.1 l2010ali 
Thomas et al.l |201 



Benson et all I2OIII: iKomatsu et al 



Riemer- S 0rensen et a 



20121 : iSanchez et alJ 



, I2OI2I: Ide Putter et all I2OI2I: IXia et all 
2012'). In this work, we use the best-ever 




measurement of the galaxy power spectrum, based on the SDSS-III 
Baryon Oscillation Spectroscopic Survey (BOSS) combined with 
other datasets of CMB, SNe and BAO to measure the neutrino 
mass, the effective number of species of the neutrinos, as well as 
other degenerate cosmological parameters taking into account the 
neutrino effects. 

We shall briefly introduce the BOSS sample in the next sec- 
tion, and describe the method adopted to model the galaxy power 
spectrum in Section 3. Before we show our results in Section 5, 
we detail the datasets included in addition to the BOSS data. We 
conclude and review the primary results of this work in the final 
section. 



2 CMASS DR9 DATA AND MEASUREMENT 



The SDSS-III dEisenstein eTai1l201lh Baryon Oscillation Spectro- 
scopic Survey J Pawson et al.l I2OI 2I) obtains spectra for two pri- 
mary galaxy samples selected from SDSS imaging data. In com- 
bination, the SDSS-I, SDSS-II, and SDSS-I II imaging surveys ob- 
tained wide-field CCD photo metrv (iGunn et a l. 1998, 2006) in five 
passbands (u, g,r,i, z; e.g.. Ipukugita et al.[ 1996h . amassing a to- 
tal footprint of 14,555 deg^, inte rnally calibrated using the 'uber- 
calibration' process described in IPadmanabhan et al] fcOOSh , and 
with a 50 percent completeness limit of point sources at r = 22.5 
(lAihara et al. 201 1). From this imaging data, BOSS has targeted 1.5 
million massive galaxies split between LOWZ and CMASS sam- 
ples over an area of 10,000 deg^ (Padmanabhan et al. in prep.). 
BOSS observations began in fall 2009, and the last spectra of tar- 
geted galaxies will be acquired in 2014. The R = 1300-3000 BOSS 



spectrographs femee et aDl2012h are fed by 1000 2" aperture opti- 
cal fibres, allowing 1000 objects to be observed in a single point- 
ing. Each galaxy observation is performed in a series of 15-minute 
exposures and integrated until a fiducial minimum signal-to-noise 
rati o, chosen to ensure a high redshift success rate, is reached 
fsee lDawson et "^1 120121 for further details). The dete rmination of 
spectr oscopic redshifts from these data is described in lBolton et alj 
(l2012h . 

We use the SDSS-III Data Release 9 fDR9:lAhn et al 



20B) 



CMASS sample o f galaxies, as analysed bylAnderson et alT 
"Reid et al.' (lOU); Samushia et al' ( 2012) : ISanchez et aP 
Nuza et al. (2012); Ross et al. (2012tlQ). It contains 264 283 mas- 
sive galaxies covering 3275 deg^ with redshifts 0.43 < z < 0.7 
(the effective redshift ^eff = 0.57. The sample is roughly split into 
two angular regions, with 2635 deg^ in the Northern Galactic Cap 
and 709 deg^ in the Southern Galactic Cap. 

We measure the galaxy p ower spectru m Pmeas(fe), using the 
sta ndard Four i er tech nique of iFeldman et al.l (1 19941) , as described 
bv lReid et all ( l2010bh . In particular, we calculate the spherically- 
averaged power in k bands of width Ak = 0.004 /iMpc~^ 
using a 2048^ grid. We use the same weighting scheme as in 
lAnderson et"aD(l2012h . lRoss et al.l(l2012bl) and Ross et al. (2012ab, 
which includes both weights to correct for the systematic relation- 
ship between target galaxy density and s tellar density and als o 
'FKP' weights, using the prescription of iFeldmanet aP (Il994l) . 
which amounts to a redshift-dependent weighting in our appli- 
ca tion. The p rocess of calculating weights is described in detail 
in IRoss et al. ( 2012b). The galaxy power spectrum was used in 
lAnderson et al. (12012 ) to extract and fit the Baryon Acoustic Os- 
cillations. 

To determine expected errors on the power spe ctrum, we use 
the 60 mock DR9 CMASS catalogues generated bv lManera et al.1 
(I2OI2I) to generate a covariance matrix C[A;i][/cj] quantifying the 
covariance of the power between the iih and jth k bins. The corre- 
lation between bins can be simply quantified using the correlation 
matrix Corr [ki ] [kj ] , 



Corr[A;i][A;j] = 



C[A)2] \kj 



(1) 



^cMh]c[k,][k,] 

Manera et al ] (l2012h used the initial conditions of a flat cosmol- 
ogy defined by Qm = 0.274, Qbh'^ = 0.0224, h = 0.70, Us = 
0.95, and ag — 0-8 (ma tchi ng the fiducial cosm ologies as- 
sumed in I White et aDl201ll and lAnderson et al.l l2oT2h . and gen- 
erated dark matter halo fields at redshift 0.55. These were pro- 
duced using a 2nd-order Lagrangian Perturbation Theory (2LPT) 
approach inspired by the Perturbation Theory Halos (PTHalos; 
IScoccimarro & Sh eth 2002). Galaxies were placed in halos using 
the halo occupation distribution determin ed from CM ASS ^0 mea- 
surements and the parameterization of Izheng et al.l tPOll . The 
DR9 angular footprint was then applied and galaxies were sam- 
pled along the radial direction such that the mean n{z) matched 
the CMASS n{z), thereby providing 600 catalogues simulating the 
observed DR9 CMASS sample. See Manera et al. (2012) for fur- 
ther details^_^ 

IRoss et alJ (l2012bl) and lRoss et"aD (l2012al) consider remain- 
ing systematic uncertainties in the power spectrum measurements, 
caused by observational effects 0. They suggest that the most con- 
servative treatment is to allow for a free parameter in the measured 



^ The systematic uncertainties only affect the power spectrum measure- 
ments on large scales (k < 0.02 h Mpc"-*^). 
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Figure 1. The power spectrum measured from the CMASS sample from the Northern (top left panel), Southern (top middle panel) and the Northem+Southem 
(top right panel) sky. The correlation coefficients obtained from the mock catalogues and the expected window functions are shown in the middle and bottom 
row. See text for more details. 



data such that 

-Pmeas(^) — 



(fc) - S[P, 



meas.nw 



(k) 



s,w(^)], (2) 



where Pmeas,w(A^) is the power spectrum measured after applying 
the weights for stellar density (the dominant known source of sys- 
tematic deviations in target density) and Pmeas,nw(A^) is the mea- 
surement without these weights. The parameter S represents an ad- 
ditional nuisance parameter that we marginalise over for all of the 
constraints we present. Note that fixing S = represents the case 
where no systematic uncertainty in the application of the weights 
for stellar density is taken into account whe n measuring the p ower 
spectrum using these weights. Following Ross et al. (2012a"), we 
put a Gaussian prior on S around zero with a stan dard deviation 0.1 , 
which is motivated by the tests using mocks in lRoss et al.l 

In the resulting measurements, data at different scales will be 
correlated as a result of the survey geometry. We must account for 
this effect when comparing theoretical power spectra, Pg(k) (cal- 
culated as described in Section 3), to Pmeas{k). We do so using 
the spherically-averaged power in the window, which we denote 
Pwin(k): the measured power spectrum Pmeas{k) is a convolution 
of the true underlying power spectrum with Pwin{k), normalised 
such that the convolved product is zero at /c = 0. For ease we con- 
vert the convolution into a matrix multiplication based on a "win- 
dow matrix" W[ki] [kj] such that 

-Pconv 

(ki) = W[h][k,]Pg{k,) - PoPwin(fc) (3) 



where 



Po = J2 W[0][kj]Pgikj)/P^i„{0). (4) 



The 2nd term in Eq. ^ is necessary because the galaxy density is 
estimated from the data itself, forcing Pmeas{k = 0) = 0. We cal- 
culate Pwin(A:) and Pg{k) in bins Ak = 0.0005 /i Mpc~\ yield- 
ing sufficient resolution to output Pconv{k) in (the measured) bin 
width Ak = 0.004 h Mdc~\ Ross et al.. C 20 12b) present a number 
of tests of this procedure based on analysis of the mock samples. 
The -Pconv (k) we determine using this approach are compared di- 
rectly t0Pmeas(A:). 

We assume a flat ACDM cosmology with Qm — 0.285, Qb — 
0.0459, h = 0.70 , Us 0.96, ag = 0.8 (approximately the best- 
fit model found bv lSanchez et al.l2()12b when calculating the power 
spectrum, and we project the scales (using Dv{z = 0.57) as the 
distance) to the cosmology to be tested. The quantity Dy is the 
average distance to the galaxies-pairs, defined as. 



DV{Z) : 



1/3 



(5) 



where Da (z) and H{z) are the physical angular diameter distance 
and the Hubble expansion rate at redshift z respectively. 

The measured power spectra Pmeas,w(A:), the correlation 
matrix Corr[A;i] [/cj] (defined in Eq. [T|) and the window matrix 
are shown in Fig. [T] (for clarity, only every 4th /c-bin 
are displayed in the middle and bottom panels). The effect of the 
survey window is clear, as the middle panels show that there is sig- 
nificant correlation between /c-bins. These correlations are larger 
in the Southern region, as is the uncertainty, due to the fact that it 
accounts for only 22% of the area observed in the total combined 
sample. 
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Parameter 


Meaning 


Prior (flat) 


Value for the vanilla model 




The physical baryon energy density 


[0.005,0.1] 


varied 


UJc = Q^ch? 


The physical dark matter energy density 


[0.01,0.99] 


varied 


Os 


100 X the ratio of the sound horizon to the angular diameter distance at decoupling 


[0.5, 10] 


varied 


T 


The optical depth 


[0.01,0.2] 


varied 


ris 


The spectral index of the primordial power spectrum 


[0.5,1.5] 


varied 


log[10iOA,] 


The amplitude of the primordial power spectrum 


[2.7,4] 


varied 


Asz 


The amplitude of SZ power spectrum when using CMB 


[0,2] 


varied 


as 


The running of the primordial power spectrum 


[ — U. ±, U. ±J 


n 

U 




The sum of the neutrino masses in the unit of eV 


[0,2] 





Neff 


The number of the neutrino species 


[1.5,10] 


3.046 


wo 


The wo parameter in the CPL parametrisation 


[-3,3] 


-1 


Wa 


The Wa parameter in the CPL parametrisation 


[-3,3] 





Qk 


The contribution of the curvature to the energy density 


[-0.1,0.1] 





r 


The tensor to scalar ratio 


[0,2] 





a 


The nuisance parameter for SN defined in Eq. (|20] 


[0.6,2.6] 


varied 


/3 


The nuisance parameter for SN defined in Eq. (|20] 


[0.9,4.6] 


varied 


bi 


The nuisance parameter for P{k) when using Eq. (|9) 


[1,3] 


varied 


b2 


The nuisance parameter for P{k) when using Eq. {9} 


[-4,4] 


varied 


N 


The nuisance parameter for P(A;) when using Eq. (|9j 


[0, 5000] 


varied 




The nuisance parameter for P(k) when using Eq. {18) 


[0.1,10] 


varied 


PUF 


The nuisance parameter for P{k) when using Eq. (fTs) 


[0, 5000] 


varied 


bQ 


The nuisance parameter for P(k) when using Eq. ([19) 


[0.1,10] 


varied 


Q 


The nuisance parameter for P{k) when using Eq. (fT9) 


[0.1,50] 


varied 


S 


The nuisance parameter for P(k) measurement systematics, defined in Eq. (|2) 


[-1,1] 


varied 



Table 1. The parameters used in our analysis and their physical meaning, ranges and values for the vanilla ACDM model. 



3 MODELLING THE GALAXY POWER SPECTRUM 

As massive neutrinos change the underlying shape of the matter 
power spectrum, resulting constraints on non-linear scales depend 
on our ability to model similar changes in shape induced by the 
physics of galaxy formation - known as galaxy bias. Various mod- 
els and parameterisations have previously been introduced to con- 
vert from the linear matter pow er spectrum to the galaxy power 
spectrum. ISwanson et alj (l2010h tried twelve different models and 
compared the resultant neutrino mass measurement using SDSS-II 
LRG sample. Models are based on two different approa ches: either 
using higher- than-linear order perturbation theory g.g. JSaito et al.l 
I2OO8 , 2009), or using fittin g formulae calibrated from A/^-body sim- 
ulations (e.g., HALOFIT: ISmith et al.|[2003l : iBird et al.ll2012h . In 
this work, we shall compare results from both approaches. 



3.1 Perturbation theory 

To extend the power spectrum to mildly nonlinear scales, a straight- 
forward approach is to include higher-order corrections to the lin- 
ear standard perturbation theory (SPT). For the system consist- 
ing of baryons and dark matter only, the higher-order corrections 
(dub bed 'one-loop' corrections) have been extensively studied 
(e.g. , Ij^s zkiewicz 198 1': 'Makino et al. 1 1 1992': 'Jain & Berts chinge3 
ll994l :1scoccimarro & Frieman 1996; Heavens et al. 1998). At the 
one-loop level of the SPT, the corrected matter power spectrum is 



plloop 
cb 



{k)^P^Uk) + P^^'\k)+P^^^\k) 



p(13)/ 



(6) 



where the subscript 'cb' denotes 'CDM plus baryons', the super- 
scripts '(22)' and '(13)' illustrate the perturbative corrections to the 
power spectrum at next- to-leading order, and the superscript 'L' 
stands for the linear matter power spectrum. With the presence of 



massive neutrinos, ISaito et al] ([2008tl2009h generalised Eq. ^ to 

PclT" (fc) = /cb^'cb°°'' (k) + (k) + flp]; (fc) , (7) 

where P^^^ [k) and Pi' (k) represent the linear power spectrum for 
total matter (CDM plus baryons plus massive neutrinos) and for 
massive neutrinos only, respectively, and the coefficient fi denotes 
the mass fraction of each species relative to the present-day energy 
density of total matter, r^m, i-e.. 



/cb 



(8) 



Q.mh'^ X 94.1eV 

A crucial assumption to derive the formula of Eq. (|7) is to treat the 
neutrino component to stay completely at linear level. ISaito et aP 
12009) shows that this assumpti on can be justified for expected 
small m ass of neutr inos (see also Shoji & Komatsul feOlO ')). Given 
Eq. Q), ISaito et al.l j2Q08, 2009) proposed a SPT-based model to 
convert the matter power spectrum to the galaxy power spectrum at 
a given scale k and a given redshift z, 

Pg{k- z) = hi [pllT^'ik- z) + 62Pb2(/c; z) + 6^Pb22(A:; z)\ +iV, 

(9) 

where bi , 62 and N denote the linear bias, nonlinear bias and the 
residual shot noise respectively, which can be derived using a SPT 
prescription (McDonald 2006). Quantities Pb2, Pb22 can be calcu- 
lated using SPT and the expression is explicitly given in Eqns. (32) 
and (36) in lSaito et aD (12009) . Note that on linear scales, Eq. ^ 
reduces to 



P,{k-z) 



blPchu{k]z) 



(10) 



Hence bi acts as a linear bias and contributes a shot noise con- 
tamina tion stemming from the stochastic bias and nonlinear clus- 
tering dHeavens et al.lll998l : [Seliak 2000: ISmith et al.ll2007h . The 
terms multiplying 62 give rise to a scale-dependent bias due to the 
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nonlinear clustering. In general, bi,b2,N vary with galaxy type, 
so we treat them as free parameters to be marginalised over in our 
analysis. 

3.1.1 Modelling the redshift space distortions 

In this section we discuss the issue of model uncertainty in the 
CMASS power spectrum, focusing on the effect of the redshift- 
space distortions (RSD). In order to quantitatively model the 
impact of the RSD on the spherically averaged CMASS power 
spectrum, we compare the following RSD models: 

RSD Model 1: Linear Kaiser At very large scales where the linear 
perturbation theory holds, the mapping from real to redshift space 
can be expressed at linear order, and the resultant redshift -space 
powe r spectrum is enhanced by the so-called Kaiser factor (iKaisej 
[l987|): 

P,^(^,/i) = (l + //i')V,(^), (11) 

where /x is the cosine between the line-of- sight direction and the 
wave vector, and / is the logarithmic growth defined as / = 
d\n.D{a) /dhia. The superscript 'S' denotes the quantity in red- 
shift space. Spherically averaging Eq. ([Til , we have the monopole 
component, 

Po(fc)= (l + Y + (12) 

RSD Model 2: Nonlinear Kaiser. The Linear Kaiser model is not 
true especially at mildly nonlinear scales (see e.g., Scoccimarro 
l2004h . To linear order, the matter density perturbation in redshift 
space can be written as, 

8^[k) ^5[k) + fiJi^e[k), (13) 

where is the divergence of the peculiar velocity field. In order to 
take into account nonlinear gravitational evolution of the velocity- 
divergence field separately, the Kaiser model is generalised to the 
Nonlinear Kaiser model as follows: 

P^{k, fi) = Pg,ssik) + 2ffi^Pg,se{k) + ffi^Pee{k), (14) 

where the galaxy density-density power spectrum, Pg^ssik) is 
modelled using Eq. (|9]), and the the galaxy d ensity-velocity power 
spectrum is modelled as jSwanson et al.ll20l"Q ), 

Pg,so{k) = bi [Pse{k) + b2Pb2,e{k)] . (15) 

Note that here we assume no velocity bias, and that the matter 
density-velocity, Pse, or the velocity-velocity Pee can be com- 
puted using perturbation theory similarly to the density-density 
one. We compare SPT with the closure a pproximation (CLA) 
(iTaruva & HiramatsulEoOSlilNishimichi et al.ll2009 ) as an example. 
The CLA is one of the improved perturb ation theories including 
the re normalized perturbation theory dCrocce & Scoccimarrol 
l2006al lbl). and the CLA power spectrum at 2-loop order is in 
an excellent agree ment with th e A^-body simulation results 
(iTaruva et al.l [20091 : ICarlson et all l2009l) . A disadvantage of the 
CLA is that it involves time-consuming integrations in the 2-loop 
order, and therefore it is c omputationally diffi cult to apply the 
CLA to MCMC analysis fsee lTaruva etal.l2012l for recent effort to 
speed up the computation). 

RSD Model 3: Nonlinear Kaiser with FoG. At smaller scales than 
the typical size of virialized clusters, the internal velocity dis- 
persion of galaxies makes the galaxy clustering pattern elongated 



along the line of sight, known as the Finger-of-God (FoG) effect 
(|jacksonlll972h . The FoG suppression is necessary for massive ha- 
los in w hich most of the CMASS galaxies exist as central galax- 
ies (e.2., iHatton & Cold [19981 with about 10% o f the CMASS 
galaxies being satellite galaxies dWhite et al.ll2oTTh . The satellites 
are expected to have larger small-scale velocity dispersion than 
central galaxies, and cause the additional FoG suppression in the 
CMASS power spectrum. In order to account for the FoG effect, 
„Scoccimarro (2004) proposed a phenomenological model in which 
the FoG suppression is described by an overall exponential factor: 

Pg{k,ii) = exp {-fayk^iJ^) 

X [Pg,ss{k) + 2ffi^Pg,se{k) + ffi^'Peeik)] , (16) 

where ay is the velocity dispersion which we treat as a free 
parameter. 

RSD Model 4: Nonlinear Kaiser with correction terms and FoG. 
Recent studies show that higher-order correlations between the 
density and the velocity divergence in the nonlinear mapping from 
real to redshift space become important to explain the redshift- 
space power spectrum especially for massive halos fra ruva etaL 
[2010; Tang et al. 201 1 ; Nishimichi & Taruva 201 1 ; Oku mura et al. 
|2012b; Reid& White 2011; O kumura et al.1 12012a). Tar uva et al. 
(12010) proposed a new model including such correction terms and 
can be generalized to biased objects: 

Pg{k,fi) = exp [-fcjYk^ii^) 

X [PgMk) + 2ffi^Pg,se{k) + fi/Pee{k) 

+bU{Kfi;p) + btB{k,fi;p)], (17) 

with /3 = f /bi. Note that terms associated with A and B include 
the linear bias dependence of bl at maximum, and bl or bf orig- 
inates from the fact that we replace / with /3 — f /bi. Also note 
that we did not include nonlinear bias terms proportional to 62 in 
the correction terms for simplicity as we expect that such terms do 
not drastically affect the discussion here. 

Now let us compare the predicted power spectra using the 
RSD models explained above, which are shown in Fig. |2] The up- 
per panel shows absolute amplitudes of the monopole power spec- 
trum (divided by the Linear Kaiser model w ith linear no-wiggle 
power spectrum from lEisenstein & Hulll998h for each model. We 
set the values of bias parameters to the best-fit values in the case 
of Am^CDM model, (61 = 2.27, 62 = 1.02, N = 2293), corre- 
sponding to the Linear Kaiser for the RSD modeling. Meanwhile, 
the lower panel plots fractional differences from the model 4 cor- 
responding to the solid black line in the upper figure. As a ref- 
erence we put the error bars taken from the diagonal components 
in the CMASS covariance matrix used in this paper. First of all, 
the Linear Kaiser using the SPT prediction (Model 1, red solid) 
overestimates the monopole amplitude at all scales compared to 
other models. The difference between the Linear (Model 1, red 
solid) and the Nonlinear Kaiser (Model 2, blue or green) is typi- 
cally - 3% at /c = 0.1 /i/Mpc or - 6% at /c = 0.2 /i/Mpc in the 
case of SPT. The Nonlinear Kaiser term does not strongly depend 
on the way Pss, Pso and Pee are computed, i.e., with SPT or CLA 
if we consider k < 0.1 h/Mpc. A comparison between the Non- 
linear Kaiser (Model 2, green) with the Nonlinear Kaiser with FoG 
(Model 3, magenta) addresses how the FoG suppress the power 
spectrum. Here we choose cry = 4.57Mpc//i which corresponds 
to the value of linear velocity dispersion. Note that, since this value 
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3.2 Fitting formula 

32.1 HALOFIT-zy 

Another approach to model the matter power spectrum in the 
weakly non linear reginie is to use fitting formulae, such as 
HALOFIT (ISmithetalJl2003h . calibrated using iV-body simula- 
tions. HALOFI T works w e ll fo r models without massive neu- 
trinos. However, iBird et alJ (Eoi2h found that HALOFIT over- 
predicts the suppression of the power due to neutrino free- 
streaming on strongly nonlinear scales, and th e discrepancy can 
reach the level of 10 percent ^Xk ^ \ h Mpc~\ lBird et al.l (l2012h 
proposed a new fitting formula, HALOFIT-z^, which is an im- 
proved version of HALOFIT when massive neutrinos are present. 
HALOFIT-i/ was calibrated using an extensive suite of A^-body 
simulations with neutrinos and it was shown to work significantly 
better than the original HALOFIT on nonlinear scales with the 
presence of neutrinos. 

To convert the matter power spectrum Phfz^ calculated using 
HA LOFIT-?^ into t he ob servable galaxy power spectrum, we fol- 
low [Swansonetal] (l20ld) combining this with free parameters for 
galaxy bias 



Pg(k\ z) = 6HFftFi/(A:; z) + Phf, 



(18) 



0.10 

k [/i/Mpc] 



where 6hf and Phf are constant model parameters representing 
the bias and the shot noise contamination respectively. 



Figure 2. An example of comparison among the RSD models. Upper, the 
monopole power spectra for the RSD models shown in the text; RSD 1: 
the Linear Kaiser {red), RSD 2: the Nonlinear Kaiser (blue for SPT and 
green for CLA), RSD 3: the Nonlinear Kaiser with the FoG prefactor {ma- 
genta), and RSD 4: the Nonlinear Kaiser plus correction terms with the 
FoG (black). Each spectrum is divided by the Linear Kaiser model with 
linear no- wiggle spectrum for clarification purpose. We consider the cos- 
mology for the CMASS mocks and the best-fit parameters of (61, 62, A^) 
in the case of AmzyCDM model. We use the linear velocity dispersion, 
cry = 4.57Mpc//i when computing the FoG prefactor. For compari- 
son, the Linear Kaiser models with 62 = —0.2 (red dashed) and with 
nil, = 0.1 eV (red dotted) are also shown. Lower: fractional difference 
of each model from the RSD model 4. The line colours and styles denote 
exactly same with those in the upper panel. We show the error bars taken 
from diagonal components in the CMASS covariance matrix as a reference. 



3.2.2 \Cole et al\tOOi) 

A n alternative to mo del the galaxy power spectrum was proposed 
bv lCole etaLl(l2005h . namely, 

PAk;z) = bl\±^Pl:Uk\z) (19) 

where k is the wave number in unit of h/Mpc and the parameters 
bq and Q change the overall amplitude and the scale-dependence 
of the linear power spectrum on small scales. It was shown that this 
fitting formula can match the redshift- space galaxy power spectrum 
derived from the simulations well if the parameters 6, Q are prop- 
erly chosen. In this analysis, we follow Swa nson et alJ (120101) to 
treat b and Q as free parameters to be marginalised over. 



is expected to be larger than that of the real CMASS catalogues, the 
FoG suppression seen in the figure could be aggressive|3 • 

In the specific case considered here, the FoG effect sup- 
presses the amplitude by ^ 5% at /c = 0.1 /i/Mpc or ^ 15% at 
k — 0.2 /i/Mpc. Finally, since the correction terms A and B in 
Eq ( FTtI essentially enhance the amplitude at scales of interest, the 
RSD Model 4 prediction (black) becomes larger than that of the 
RSD Model 3 (magenta) by - 4% at A; = 0.1 /i/Mpc or - 10% 
at /c = 0.2/i/Mpc. 



^ This argument comes from the fact that almost all of the CMASS galax- 
ies are central galaxies in the centre of massive halos and their veloc ity dis- 
persions are expected to be small (g.g. jNishimichi &' Taruvall201lh . Also, 
the effect of satellites on the monopole clustering is shown to be negligibly 
small, which is not true for the quadrupole (.Reid et al.,2012: Okumura et al.. 
[20r2ah 



4 OTHER DATASETS AND MCMC ANALYSIS 

To tighten the measurement of the summed neutrino mass, we com- 
bine the CMASS data with that from other cosmological surveys, 
including CMB, SN and Hubble parameter determinations. 

4.1 CMB and SN 

We use likelihoods based on the WMAP 7 -year CMB data set in- 
cluding the temperature-temperature and temperature-polarisation 
angular power spectra (Larson et al. 201 1 ). 

There are two recently compiled, and publicly available sets 
of SN that are co mmonly used to se t cosmological contraints - 
the SNLS 3-year (Iconlev et all I20T1I) and the Union2.1 sample 
dSuzuk i et al.' 2012). We chose to use SNLS3 data in this work be- 
cause it includes a more homogeneous sample of SN at higher red- 
shift, which may make it less susceptible to systematics errors. We 
have also performed our analysis for some cases (e.g., weighing 
neutrino mass assuming a ACDM background cosmology) using 
Union2.1 data for a comparison (see Fig. |6j. For the SNLS data. 
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z 




A{z) 


6dF 


0.106 


0.336 ±0.015 




SDSS-TT 


0.20 
0.35 


0.1905 ±0.0061 
0.1097 ±0.0036 






0.44 




0. 17 1 ± 0.034 


WiggleZ 


0.60 




0.442 ± 0.020 




0.73 




0.424 ±0.021 


BOSS 


0.57 


0.07315 ±0.00118 





Table 2. The BAO measurements used in this work. The quantity Vs/Dy 
is the ratio of the sound horizon io Dy defi ned in Eq Js}' ^i^d A{z) is 
the parameter defined in lEisenstein et al.l (l2005b . 



the measured ap parent magnitude mmod is dConlev et alj |201 ll : 
iRuiz et alJlioTi). 



mmod - 51ogio ( -^<^L 



■a{s-l)+pC + M, 



(20) 



where di, is the luminosity distance, and the presence of nuisance 
parameters a and /3 allows the inclusion of the possible degeneracy 
between mmod, 'stretch' s and color C. This is necessary because 
in principle the peak apparent magnitude of a SN la is degenerate 
with the broadness (stretch) of the light curve (b roader is brighter) , 
and with the colour as well (bluer is brighter) (iRuiz et alj|2012l) . 
Thus in the likelihood calculation, where we essentially contrast 
the theoretical prediction to the measured mmod for each SNe and 
calculate the usual quadratic sum of the signal-to-noise ratio, we 
numerically marginalise over a and {3, and analytically marginalise 
over the constant magnitude offset M. 



4.2 BAO 

We also combine th e baryonic acoustic oscillation s (BAO) mea- 
surem ents from 6dF feeutler et al."20 1l|), SDSS-I I dPercival et alJ 
|201Q|) and the WiggleZ survey (Blake et alJboilh . The data used 
is shown in Table d Since the mea surements of CMASS P(k) and 
BOSS BAO ( And erson et all 12012) use the same galaxy sample, 
we never use the P{k) and BAO measurements of BOSS simulta- 
neously. But we shall compare the neutrino mass constraint using 
BOSS P{k) and BAO respectively (see Fig. [3] and later discus- 
sions). The difference in the constraints will quantify the degree to 
which it is the growth supression effect of massive neutrinos that 
provides the information we recover. 



4.3 Parameterising the Universe and the MCMC engine 

The set of cosmological parameters that we will simultaneously 
measure is given in Table [T] which also provides the basic defi- 
nitions, the prior assumed, and the values for the vanilla ACDM 
model for these parameters. We use Gs rather than Hq to parame- 
terise the background expansion o f the universe since is less 
degenerate with other parameters (iLewis & Bridlell2002b . A few 
require lengthier definitions, including wq and Wa, which param- 
eterise the equation-of- state (EoS) parameter w of dark energy as a 
function of the scale factor a 



This is the so-cal le d CP L pa rametris a tion proposed by 
IChevallier & Polarskil (l200lh and iLinded (l2003h . Parameters 
Ag,ng and as p arameterise the primordial power spectrum Pin 
jKosowskv & Turner 1995 ), 

lnPin(A;) = \nAs + [ns(k 



2 



In 



(22) 

where /cpv denotes the pivot scale for the parametrisation and we 
choose /cpv = 0.05 h Mpc~^ in this work, 



We use a modified version of CAMB (iLewis et alJl200Qh to 
calculate the observables and compute the likelihood by compar- 
ing to observations. In our general parameter space (see Table [TJ, 
w{a) evolves with the scale factor a, and it is allowed to cross w — 
— 1. T his means that we allow the quintessence (Peebles & Ratral 
ll988l:[Ratra & Peebles 1988), p hantom ( Caldwell 2002), k-essence 
(lArmendariz-Picon et alj I'200Q|) and quintom dPeng et all l2005h 
dark energy models in the global fitting. Note that when w{a) 
crosses the divide, i.e.. 



wq < —1 and wq + Wa > — 1, 
or wq > —1 and wq -\- Wa < — 1, 



(23) 



ri(a) = Wo + Wa{l — a). 



(21) 



dark ener gy fluid has multi-components dPeng et al.l l2005l: 
IhuI I2OO5I) . In this c ase, w e follow the prescription pro- 
posed by IZhao et aP (l2005h to calculate the evolution of 
dark energy perturbations [j. We use a modified version of 
COSMOMC (^ewis & Bridle 2002), which is a Monte Carlo 
Markov Chain (MCMC) engine, to efficiently explore the multi- 
dimension al parameter space using the Metr opolis-Hastings algo- 
rithm (Metr opolis eTaDll953l:lHastingslll970h . 

To summarise our default assumptions, unless specifically 
mentioned, 

• we use Eq. ^ to model the galaxy power spectrum; 

• we use WMAP seven-year power spectra and the SNLS three- 
year measurement for the CMB and SN data respectively; 

• we use the CMASS galaxy power spectrum measurement 
truncated at /c = 0.1 /i Mpc"^ for the 'CMASS' data; 

• when using the 'CMASS' data, we always combine the BAO 
measurement shown in TdbXo^except for the BOSS measurement; 

• all the neutrino mass constraints are the 95 percent CL upper 
limit; 

• all P{k) shown are evaluated at the mean redshift of BOSS, 
which is Zeff = 0.57. 



5 RESULTS 

5.1 Neutrino mass measurement in a ACDM universe 

We now measure the summed neutrino mass T^rriu assuming that 
the background cosmology follows that of a ACDM model, i.e., 
only the vanilla cosmological parameters cj^, cjc, Os, r, ris. As, 
T^rriu, and the nuisance parameters are allowed to vary. 



5.1.1 The choice of /cmax 

Fitting to a wider range of scales has the potential to provide 
better constraints on T^rriu, as we have a longer lever-arm with 



^ Note that alternative approaches hav e been proposed and can yield con- 
sistent result (e.g., see lFang et al.ll2008h . 
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Figure 3. The purple and green contours on the top layer in left panels: the 68 and 95 percent CL contour plots for neutrino mass and VLyi obtained from the 
joint dataset including CMB+SN+CMASS power spectra cut off at various k illustrated in the figure; The blue contours on the bottom layer in left panels: 
the 68 and 95 percent CL contour plots for neutrino mass and Vlyi obtained from the joint dataset including CMB+SN+CMASS BAO; Right panel: the 
corresponding ID posterior distribution of neutrino mass. A ACDM model is assumed for the background cosmology. 
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Figure 4. The CM ASS data used in the analysis and the best fit power 
spectrum assuming a AmiyCDM cosmology. The data and spectra are both 
rescaled using the linear matter spectrum for the best fit model. The upper 
and lower panels show the cases of A^max = 0.1 and 0.2 /i Mpc"-*^ respec- 
tively. 



bias of the CMASS galaxies means there is little information over 
the range 0.1 < k < 0.2 /iMpc~^, assuming our model is appro- 
priate on these scales is not worth the added risk, and we choose 
A:max = 0.1 /iMpc~^ as our default. All of the results presented 
are based on this conservative limit. 

Fig. |4] shows the goodness-of-fit by over-plotting the best-fit 
model on top of the observational data for the cases of /cmax = 
0.1, 0.2 /iMpc~^. The quantity shown is the ratio between the 
galaxy power spectrum and the linear matter power spectrum, 
Pg {k)/P^i,^ (see Eqns. Q and (Igjl). From Eq. ([TOl ), we can see that 
on large scales, this ratio is roughly the linear galaxy bias squared 
bf. Our MC MC analys is suggests bi 2, which is consistent with 
the result in lReid et al] (12012) . 



5.1.2 The choice of the galaxy modelling 

Next we shall test the effect of the choice of the galaxy modelling. 
To test the limiting scales to be fitted in the last section, we used 
Eq. ^ to model the galaxy power spectrum. In Fig. |5] we show 
the contours for Y^rr iu and for the HALOFIT-zy (Eq.[T8]) and 
the lCole et al.l (l2005h (Eq.[T9l) models. The difference between con- 
straints calculated assuming these three models is marginal, which 
is reasonable as they only differ in form for /cmax > 0.1 Mpc~^. 



which to assess the contribution from the neutrinos. However fit- 
ting to smaller scales relies more heavily on the model for galaxy 
formation. In Fig. [5] we compare results fitting to scales with 
A:max = 0.1, 0.2 /i Mpc~^ when all data are combined (purple and 
green contours). As shown, the change is marginal, with the up- 
per limit for llmi, only lowered to 0.338 eV from 0.340 eV when 
A^max is increased from 0.1 to 0.2 h Mpc~^. This is understandable 
given our default galaxy bias model: when /cmax is larger, the non- 
linear P{k) data simply constrains the nuisance parameters 62 in 
Eq. ^ rather than T^rriu. We should expect that, for larger /cmax, 
the model Eq. ^ becomes less reliable since it is based on pertur- 
bation theory. Given that our lack of knowledge of the non-linear 



5.1.3 The choice of SN data 

The left panel of Fig.[6]shows 95 percent CL contour plots for T^rrit^ 
and comparing various data combinations. Comparison with 
the left panel shows the effect of including the Union2.1 rather 
than SNLS3 SN data. SNLS3 data provides tighter constraint on 
Qm, although the measurements of the neutrino mass are similar, 
namely, 

Sm^ < 0.340 eV (WMAP7 + SNLS3 + CMASS), 
Em,. < 0.334 eV (WMAP7 + Union2.1 + CMASS).(24) 

In both cases, we see that CMASS data help to reduce the allowed 
parameter space dramatically. 
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Figure 5. The 68 and 95 percent CL contour plot for neutrino mass and Qm using three different galaxy modelling. 
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Figure 6. The 95 percent CL allowed parameter space for (Qm^^^tiu) 
obtained from the joint dataset including CMB (red), CMB+SN(purple), 
and CMB+SN+CMASS P{k) (white) in a Amz^CDM model. The SN 
data of SNLS 3-year and the Union2.1 sample were used for the left and 
right panels respectively. The CMASS P(k) was truncated at /cmax = 
0.1/iMpc-^ 

5.1.4 The effect of the redshift space distortion modelling 

We now test how the neutrino mass measurements are affected by 
the choice of the RSD model. In Figs [T] [8] and Table [S] we show 
the constraint on neutrino mass (and related parameters) derived 
from the same dataset (the CMASS P{k) is cut at /c = 0.1/i/Mpc) 
using different RSD modelling. As we can see, the neutrino mass 
constraints are generally similar. This is mainly because the effect 
of scale-dependent RSD up to A; = 0.1/i/Mpc is rather mild and 
is unlikely to be degenerate with the neutrino mass. Notably, we 
do not find a significant degeneracy between the neutrino mass and 
the FoG parameter, even though they both suppress the amplitude 
of the power spectrum. Given these results, in this work we adopt 
the RSDl model, which is equivalent to the SPT model (Eq.[9]), as 
a default just for simplicity. 

5.7.5 P{k)orBAO? 

A strong feature of the galaxy power spectrum is the BAO signal, 
whose location depends on Vs/Dy, where Vs is the sound horizon 
at the baryon drag epoch, and Dv is the average distance to the 
galaxies pairs defined in Eq (|5]). The full CMASS P{k) includes 



T 1 1 1 r 




2 m reV] 



Figure 7. The 1-D posterior distribution of the neutrino mass using different 
RSD modelling as illustrated in the legend. 







[eV] 




RSD modelling 


95 percent CL 


Mean 


(7 


RSDl 


< 0.340 


0.149 


0.104 


RSD2 


< 0.336 


0.140 


0.103 


RSD3 


< 0.384 


0.158 


0.118 


RSD4 


< 0.324 


0.135 


0.099 



Table 3. The neutrino mass constraint for various RSD modelling. 

information on additional physical processes, including those re- 
lated to neutrinos. In some cases, e.g., for dark energy parameters, 
BAO and the full P (k) can provide similar constraints (see, e.g., 
ISanchez et al.ll2012h . By comparing our constraints using the full 
P{k) to the BAO only result, we can determine the amount of in- 
formation about physical processes related to neutrinos that is en- 
coded in the CMASS P{k). 

For the constraints on the summed neutrino mass. Fig. |3] and 
TablelH show that the constraint from including the CMASS BAO 
measurement, rather than the full power spectrum is much weaker. 
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Figure 8. The 1-D and 2-D constraint on neutrino mass and related parameters. 



namely, 

Sm^ < 0.579 eV (WMAP7 + SNLS3 + CMASS BAO), 
Em, < 0.340 eV (WMAP7 + SNLS3 + CMASS P{k)). 

Clearly, the broadband shape of the CMASS P{k) contains signif- 
icant information on physical processes related to neutrinos, par- 
ticularly on the small-scale damping of P(k) caused by the free- 
streaming of neutrinos. 

5.1.6 Summary 

The measurements of the neutrino mass in a ACDM background 
cosmology are summarised in Table |4] We find that 

(i) As long as WMAP7 and CMASS P{k) data are used, us- 
ing SNLS3 or Union2.1 SN data yield similar neutrino mass con- 
straints; 

(ii) The neutrino mass constraint cannot be improved when fit- 
ting to scales /cmax ^0.1 /iMpc~^; 

(iii) If /cmax = 0.1/iMpc~^, three galaxy models for the 
galaxy power spectrum described in Eqns. (|9), dlS] and ([19), and 
the 4 models for redshift space distortions described in section [TT] 
yield similar results; 

(iv) The constraint is significantly degraded if the BAO only are 
used instead of the full power spectrum. 

5.2 Neutrino mass measurements in more general 
background cosmologies 

In this section, we reinvestigate the summed neutrino mass mea- 
surement in more general cosmologies, when several parameters in 
Table [T] are allowed to vary. Figs. |9] and (Tolpresent likelihood con- 
tours showing neutrino mass constraints in the models considered 

5.2. 1 Neutrino mass in a wrriu CDM universe 

Neutrino mass and w (assumed to be a constant here) of dark en- 
ergy are generally degenerate because both can affect the shape of 



Floating parameters 


Acronym 


V 


ACDM 


V + T^rriu 


Ami. CDM 


V + Em^ + A^eff 


AA^efF^i/CDM 


V + Emi/ + WQ 




V + Em^, +WQ+Wa 


WQWaTTli^CDM 


V + Hrriiy + OLs 


Ami, as CDM 


V + WQ 


k;CDM 


V + Wq + Wa 




V + VLk 


oCDM 




A^sCDM 


V + r 


ArCDM 


V + Hm^ + A^eff +WQ+Wa+^K + r + OLs 


All float 



Table 5. The floating cosmological parameters and the corresponding 
acronym for each models. The symbol V stands for the set of default vanilla 
parameters defined in Table [T] 



the matter and CMB power spectra (iHannestadl l2005h . We indeed 
see this degeneracy in Fig. [TT] and this is why the neutrino mass 
constraint is diluted if w is allowed to float, as shown in panel (A) 
in Fig. [21 In Table [6] we can see that when w is varied, the 95 
percent CL upper limit for neutrino mass is relaxed from 0.340 to 
0.432 eV, a 27 percent degradation. 



5.2.2 Neutrino mass in a woWaTriuCDM universe 

The dynamics of dark energy can further dilute the neutrino mass 
constraint (Xiaetal. 2007). We measure the neutrino mass again 
by parameterising dark energy using the CPL parametrisation, 
Eq. ([2D. As shown in panel (B) of Figs. [9] and [TO] and Table [6] 
the neutrino mass constraint is further diluted to 0.465 eV, which is 
a 37 percent degradation. 
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[eV] 




Galaxy modelling 


Data 




95 percent CL 


Mean 


cr 


SPT, Eq. H 


CMASS P(fc),/Cmax = 


0.1/iMpc-i 


< 0.340 


0.149 


0.104 


SPT, Eq. 12} 


CMASS P(fc),A^max = 




< 0.338 


0.152 


0.097 


HALOFIT,Eq. ([18} 


CMASS P(/e),fcmax = 


0.1/iMpc-i 


< 0.388 


0.167 


0.117 


HALOFIT,Eq. {TUl 


CMASS P(/C),fcmax = 


0.2/iMpc-i 


< 0.402 


0.173 


0.121 


Cole etal. (2005). Ea. (fT9l 


CMASS P(/C),fcnaax = 


0.1/iMpc-i 


< 0.399 


0.168 


0.121 




CMASS BAG 


< 0.579 


0.261 


0.172 



Table 4. The neutrino mass constraint for various galaxy modelling and data choices. 




Figure 9. Green contours in the back layer: the 68 and 95 percent CL con- 
straints for the summed neutrino mass and flyi assuming different cosmolo- 
gies shown in the legend. Transparent magenta contours on the front layer: 
the constraint assuming a Amz^CDM cosmology for a comparison. 

5. 2. 3 Neutrino mass in a ANef^rriu CDM universe 

The effective number of the relativistic species A^eff is directly 
related to Q^mh^ , the physical matter density, and Zeq, the red- 
shift of the matter-radiation equality. Consequently, varying A^eff 
can change the ratio of the first to the third peak of the CMB 
temperature-tem perature spectrum, an d also the shape of matter 
power sp ectrum dKomatsu et alj|201lb . The standard value of A/eff 
is 3.046 (iMangano et alJllOOSb . but in general, A^eff can be consid- 
ered a free parameter to be constrained by data. The bound of A^efF 
can be used to study the neutrino physics, and investigate the upper 
bound of energy density in primordial gravitational waves with fre- 
quencies > 10^^ Hz dKomatsu et alj|201lh . As we see from panel 
(C) in Fig. [9] and in Table[6l relaxing A'eff dilutes the neutrino mass 
more than the dark energy parameters: the constraint is relaxed to 
0.491 eV, which is a 44 percent degradation. 



' 1 ' 1 
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wm CDM 
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+CMASS 
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EmJeV] 

Figure 11. The 68 and 95 percent CL contour plots for w and neu- 
trino mass obtained from the joint dataset including CMB-i-SN (green) and 
CMB-i-SN-i-CMASS power spectrum (magenta) in a wrrii^CDM cosmology. 



5.2.4 Neutrino mass in the most general cosmology 

Finally, let us consider the most conservative case in which all the 
cosmological parameters listed in Table [T] are allowed to vary si- 
multaneously: in addition to wq, Wa, & A'eff, the curvature Qk, 
running as and tensor- to- scalar ratio r are varied. All of these pa- 
rameters can, in principle, be degenerate with the neutrino mass. 



e.g., as can mimic the effect o f neutrinos on P (k) by titling the 
power, thus changing the shape dFeng et al.l2006r) . As expected, the 
neutrino mass constraint is weakened significantly to 0.879 eV, a 
160 percent degradation compared to that in a ACDM background 
cosmology. 
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ACDM 


Amz^CDM 


Amz^A^effCDM 


wniiyCDM 


WoWaTTluCDM 


AniuasCDM 


All float 


T^ruu [eV] 





< 0.340 


< 0.491 


< 0.432 


< 0.618 


< 0.395 


< 0.821 


A^eff 

WQ 

Wa 


3.046 

-1 



3.046 
-1 



4.308 ±0.794 
-1 



3.046 
-1.081 ±0.075 



3.046 

-0.982ig:l^^ 

-U.a«_Q ggg 


3.046 

-1 



4 ns2+0-8'^o 
-0 964+°-^'^^ 

PI 701 +0.982 
^- ''^-'^-1.053 






















-0.411 ± 1.02 


r 




















< 0.440 


OLs 

















-0.012 ±0.019 


-0.032 ±0.030 


lOOl^c/i^ 


2.265 ± 0.050 

n SI 7+0-243 


2.278 ± 0.054 
11 217+0-254 

-L-L-^J- ( -0.256 


2.261 ±0.051 

1 S SI +1-32 


9 9«Q+0-050 
z.zuo_Q Q52 
11 455+0.352 


9 9/1-1+0.053 
z.z^j__Q Q52 

11.765 ±0.472 


2.246 ± 0.060 
11 sos+O-264 


9 9fiQ+0.080 

Z.ZUO_Q Qg]^ 

12 974+1-304 

-^^•^'^-1.329 


r 


1 nSQ7^ 
i.UCiy /_Q QQ25 

0.0870lg;°°?g 


1 nsQ4^ 

J_.UO»^_Q QQ27 

^ ^^^^-i-D nnfi5 


1.0350 ±0.0038 
0-088lio:oo69 


1.0390 ±0.0026 

U.U004_Q QQg3 


1.0387 ±0.0026 

^ ^^--^-1-0 0061 

0.0850lO;00^^ 


1.0395 ± 0.0026 

^ -1-0 0071 

U.UyZ±_Q QQgQ 


1.0364 ±0.0039 

^ ^^^^^-1-0 006Q 


Tls 


0.967 ±0.012 


0.969 ±0.012 


0.982 ±0.014 


0.966 ±0.012 


0.961 ±0.013 


0.954 ±0.026 


0.965+°°« 

— U.U4Z 




3.085 ±0.033 


3.076in?^ 


3.115 ±0.038 


3.078 ± 0.033 


3.081 ±0.033 


3.084 ± 0.036 


3.115 ±0.042 


lOO/i 


0.2796 ± 0.0097 


RQ 40+1-03 

Oy.4U_Q gg 


n 9S9Q+0-0100 

74 14+3-17 
'^•-^^-3.14 


n O7f?7+0.0108 

70.4611:^3 


n 2815+°-°^^° 

U.Z5iO_Q Q-L28 

70.56 ± 1.43 


69.33ll;°5 


0.2798t°o.Slf4 
73.78+11^ 


erg 


0.816 ±0.021 


--+0.036 
^- ' ' ^-0.034 


0.807 ± 0.044 


0.786 ±0.046 


0.792 ± 0.053 


0.772 ± 0.039 


796+0-063 

U. 'yD_o.064 


Age [Gyrs] 


13.79 ±0.10 


13.84 ±0.12 


1 2 Q'^+0-53 


13.86 ±0.12 


13.84 ±0.12 


13.83 ±0.11 


13.28 ±0.59 



Table 6. The mean and the 68 percent CL errors for the cosmological parameters measured assuming various models. For the parameters with a one-tailed 
distribution, such as the neutrino mass and r, the 95 percent CL upper limit is presented instead. 



Minimal model 
Parametrisation Constraint [CMASS P{k)] 



Parametrisation 



All float 
Constraint [CMASS P(k)] 



Constraint [CMASS BAG] 



^rrijy [eV] 


Ami^CDM 




AA^eff^i^CDM 




WQWaCDM 


Wa 


WQWaCDM 


IOOQk 


oCDM 


r 


ArCDM 


as 


Aa^CDM 




ACDM 




ACDM 


Ss 


ACDM 


T 


ACDM 


Us 


ACDM 


ln[10iOA,] 


ACDM 




ACDM 


lOO/i 


ACDM 


Age [Gyrs] 


ACDM 



< 0.340 
4.308 ± 0.794 

-1.041 ±0.143 
-0.111 ±0.708 

-0 264+0-466 

< 0.198 
-0 017+°-°^^ 

-0.017 



AU float 
AU float 
AU float 
AU float 
AU float 
AU float 
AU float 



< 0.821 

4 ns2+0-8'^o 

4.U.5Z_Q gg4 

-0 964+°-^'^^ 
u.yo^_Q j^gg 

n 701 +0.982 
-U. '^<3i_i.o53 

-0.411 ± 1.02 

< 0.440 
-0.032 ±0.030 



< 1.143 
4.324 ±0.881 
-0.899t°jf, 

U.DOi_Q ggg 

< 0.405 
-0.023 ± 0.030 



2.265 ± 0.050 

-I -I 09+0.243 

ii.OZ_Q 245 

i.UOy^ _Q QQ25 

n n«7n+0-0063 
U.U87U_o.oo70 
0.967 ±0.012 
3.085 ± 0.033 



AU float 
AU float 
AU float 
AU float 
AU float 
AU float 



2.268Z 



n + 0.080 

^-0.081 

12 Q74+1-304 
iZ.y 329 

1.0364 ±0.0039 
0927+0-0069 

u.uyz I _o.o082 

965+0-043 
u.yoo_Q Q42 

3.115 ±0.042 



cy .+0.082 

Z.ZO^_Q Qg3 

1 S 671 +1-587 
i.3.D^i_;L.575 

1.0351 ±0.0038 
0922+°-°°^^ 

U.U»ZZ_Q QQ-^g 

971+0-042 
u.y^ -L_o.041 

3.121 ±0.045 



0.2796 ± 0.0097 

69.72ig:^g; 

13.79 ±0.10 



AU float 
AU float 
AU float 



0.2798; 



73.78 



+ 0.0146 
-0.0144 
+3.16 



-3.17 



13.28 ±0.59 



74.561^3:^4 
13.20lg:^g^ 



Table 7. The constraint of the cosmological parameters in the optimistic (vanilla parameters with a minimal extension to include the parameter concerned) 
and the conservative case where aU parameters are allowed to float. For the conservative case, we show the result using the CMASS P{k) and CMASS BAG 
respectively. 



5.3 Measurement of other cosmological parameters 

In this section, we present the measurements of other cosmological 
parameters. For each parameter x, we first present the constraint 
in the xCDM modeQ We then provide the constraint in either the 
xmzyCDM, or the 'All float' models, according to the degeneracy 
between x and the neutrino mass. In this way, we can see the im- 



4 If x is a member of the set of the vanilla parameters presented in Table[T] 
we present the constraint in ACDM. 



pact of neutrino mass, and other cosmological parameters on the 
measurement of parameter x. 

5.3.1 Qm and Hq 

In ACDM, Qm and Ho are well determined, namely, 

Om = 0.2796 ± 0.0097, Ho = 69.7210.9? km/s/Mpc (25) 
If the neutrino mass is marginalised over, the constraint becomes, 
^^M = 0.2804lg:glg^, ifo = 69.40io.99 km/s/Mpc (26) 
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)5 I ' ^ ' ^ ' ^ ' ^ ' ^ ' ^ 

0.20 0.22 0.24 0.26 0.28 0.30 0.32 



^^M Figure 14. The 68 and 95 percent CL contour plots for w and ob- 

tained from the joint dataset including CMB (blue), CMB+SN (green), 
CMB+SN+CMASS power spectrum (magenta). The neutrino mass is fixed 
Figure 12. The 68 and 95 percent CL contour plots for Qm and if o ob- to zero and allowed to vary in the left and right panels respectively, 

tained from the joint dataset for A CDM (black solid), Aruiy CDM (ma- 
genta filled) and 'AH float' models (green filled). 



0.04 




Figure 13. Left: green contours on the back layer: the 68 and 95 percent 
CL constraints for fl^ and Qm in the All float' cosmology. Transparent 
magenta contours on the top layer: the same for the oCDM cosmology; 
Right: the corresponding ID posterior distribution of fix- The vertical 
black dashed line shows Qk = to guide eyes. 

If all parameters are marginalised over, the constraint is diluted to, 

Qm = 0.2798io.gl36, Ho = TS.TSla.lr km/s/Mpc (27) 

From these numbers and also Fig. [12] we can see that marginalis- 
ing over neutrino mass can only dilute the constraint on these back- 
ground parameters by roughly 10 percent, but if other parameters 
including wo^Wa, vary, the degradation is significant. 

5.3.2 Curvature Qk 

The {^k,^m) contour and the posterior distribution of Qk 
are shown in Fig. [13] In the 'All float' model, the Qk con- 
straint is weakened by more than a factor of 2, namely, Qk = 
-0.002641° °°^^^ (oCDM),andllK = -0.00411 ±0.01029 (All 
float). But in both cases, Qk is very consistent with zero. 

5.3.3 The constant w 

Due to the degeneracy between w and the summed neutrino mass 
shown in Fig. \TT\ the constraint of w is slightly relaxed when 



1 


quintom ^^^^s^ 


^■AU float 

1 WgW^mJODM 







W quintessence - 


-1 


- ^ 




phantom \\ 




-2 




v^Wy \\ quintom 


-3 


■ 





-2 -1 



W 

Figure 15. The 68 and 95 percent CL contour plots for wq and Wa obtained 
from the joint dataset in three cosmologies: WQWaCDM (black curves), 
WQWaTniyCDM (magenta) and 'AH float' (blue). The black dashed lines 
subdivide the parameter space according to the theoretical dark energy 
model predictions, and the yellow star illustrates the ACDM model. 

the summed neutrino mass is marginalised over, namely, w = 
-1.068 ± 0.072 (^CDM); w = -1.081 ± 0.075 (^m^CDM). 
If all the parameters are marginalised over, the constraint is further 
weakened to ^ = -1.067 ± 0.121. 

5.3.4 Wo and Wa 

Similarly, the CPL parameters wo and Wa are less constrained when 
the other cosmological parameters, including the neutrino mass, are 
simultaneously fitted: 

^0 = -1.041 ± 0.143, Wa = -0.111 ± 0.708 (woWaCDM), 
wo = -0.984 ± 0.157, Wa = -0.704 ± 0.901 (wo^am^CDM), 
^0 = -0.964 ± 0.168, Wa = -0.731 ± 0.970 (All float). 

The contour plot for wo , Wa is shown in Fig.[T5]for three cosmolo- 
gies. As we can see, marginalising over the neutrino mass can sig- 
nificantly dilute the measurement of wo and Wa, and also shift the 
contour towards a more negative Wa • This is understandable - when 
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0.15 
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^ — All float 
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Figure 16. Left panel: the 68 and 95 percent CL contour plots for A^eff 
and flyih^ obtained from the joint dataset for AA^eff^i^CDM (magenta) 
and 'AH float' models (green); Right panel: the corresponding ID posterior 
distribution of A^gff • 



the massive neutrino, whose equation-of- state is non-negative, is 
assumed to exist in the universe, a dark energy component with 
a more negative equation-of- state (assumed constant) is needed to 
compensate. This is why we see an anti-correlation between w and 
the neutrino mass in Fig[TT] If the dynamics of dark energy, repre- 
sented by a nonzero Wa for example, is allowed, the neutrino mass 
anti-correlates with the effective equ ation-of- state Wef^ of dark en- 
ergy (see Eq(57) in lzhao et al.ll2005h . As a result, either wo or Wa 
needs to be more negative if 12m u / 0. 

5.3.5 Nei, 



The number of neutrino species is measured to be 4.308 ± 0.794 
and A^eff = 4.412ig;^^^, in the AA^effm^CDM and 'All float' 
models respectively. In both cases, the standard value A^eff = 3.046 
is only consistent with this at the 1.6a level. 

Fig. [16] shows the contour plot of A^eff and Qmh^, and the 
posterior distribution of A^eff for both ANef^rriuCDM and 'All 
float' models. As we see, A^eff strongly correlates with Qyih^ as 
expected since A^eff is directly related to Qyih^ (see Eq. (53) in 
iKomatsu et aDl201 ll) . 



5.3.6 Tensor mode and the inflation models 

Various inflation models can be classified on the {ris , r) plane. 
For example, the Small-field and Large-field models predict r < 
1(1 - Us), Us ^ 1 and |(1 - Us) < r < 8(1 - n^), Us ^ 1, 
respecti vely, while Hybrid inflation predicts r > |(1 — ris) and 
Us free dPeiris et aLlbOOSl i pTegmark et al.ll2004b . To test various in- 
flation models, we show the {ris , r) contour plot in Fig. [T7] The 
constraint is 



= 0.982 zb 0.012, r < 0.198 (ArCDM), 
r < 0.440 (All float). 



ris = 0.965_o.o42 



(28) 



We also show the prediction of two Large-field inflation models, 

, ^ _ 8 ^ _ 1 2 



Us 



N ' 

= 1 - ^ 



where 0, V denote the inflaton and the potential of the inflaton re- 
spective, and N stands for the e-fold number. We show the model 
predictions for both the 0^ and 0^ cases for A^ = 50, 60. As we can 
see, in the ArCDM cosmology, the 0^ model is significantly dis- 
favoured, but the (j)^ model is still allowed. However, both models 
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Figure 17. 68 and 95 percent CL contours for (us , r) in the ArCDM (ma- 
genta) and the 'AH float' models (green). The solid black lines subdivide 
the parameter space to show different inflation model predictions. The blue 
dashed and red dash-dotted lines illustrate the models of (jP' and 0^ models 
respectively. The small and large dots on each line show the model predic- 
tion for the e-fold A^ = 50 and 60 respectively. See text for more details. 



0.06 
0.04 
0.02 
0.00 
-0.02 
-0.04 
-0.06 
-O.OJ 




1.85 0.90 0.95 1.00 1.05 
n 



1.10 



1.0 
0.8 
0.6 
0.4 
0.2 
O.Q 




0.10 



Figure 18. Left panel: 68 and 95 percent CL contours for Us and as ob- 
tained from the joint dataset for Kcxs CDM (black solid), Aasrriiy CDM 
(magenta filled) and 'AH float' models (green filled); Right panel: the cor- 
responding ID posterior distribution of as. 



are consistent with observations for the case where all the param- 
eters are allowed to vary. The scale-invariant spectrum (us = 1) 
is outside the 68 percent CL region of the resulting contour, but is 
within the 95 percent CL region in both cases. 



5. 3. 7 Running of the primordial power spectrum 

The running of the primordial power spectrum, as, is constrained 
to be as = -0.017lo:oi7 in the Aa^CDM cosmology. S ince both 
as an d the neutrino mass can change the shape of P{k) CFeng et al.l 
l2006h . we redo the analysis with the neutrino mass also allowed to 
vary. The resultant measurement is slightly diluted, namely, as = 
—0.012 zb 0.019. When all the parameters vary, the constraint is 
further weakened to be as = —0.032 ± 0.030. The contour plot 
for Us and as and the ID posterior distribution of as are presented 
in Fig. [18] As shown, as is consistent with zero in all of the cases 
we consider. 
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Figure 19. The ID posterior distribution of the cosmological parameters Hsted in Tableland the 2D 68 and 95 percent CL contours within them. The shaded 
contours and blue solid curves show the results calculated using the CMASS power spectrum, while the unfilled contours and the red dashed curves show the 
result using the CMASS BAO distance measurements. The CMB, SN and other BAO datasets are combined for both cases. The black vertical lines show the 
standard values for some of the parameters. 



^rrijy [eV] 



Reference 



Galaxy data used 



< 0.80 (0.81) 

< 0.51 

< 0.26 (0.36) 

< 0.29 (0.41) 

< 0.48 (0.63) 

< 0.32 

< 0.340 (0.821) 



Saito et al. (2009) 
Sanchez e t al. (2012) 
de Putter et al. (2012) 
Xia et al. (2012) 
Wang etTl. (2012) 
Parkinson et al. (2012) 
This work 



3D power spectrum of SDSS-II LRG 
SDSS-III CMASS two-point correlation function combined with BAO of other surveys 
SDSS-III DR8 LRG angular power spectrum (Aihara et al. 2011; Ho et al. 2012) 
Angular power spectrum of CFHTLS galaxy counts 
Weak lensing measurement of the CFHTLS -T0003 sample 
Angular power spectrum of WiggleZ galaxy counts 
SDSS-III CMASS 3D power spectrum combined with BAO of other surveys 



Table 8. Recently published neutrino mass constraints found in the literature. For the neutrino mass constraints, the number in the parentheses shows the 
measurement obtained under more conservative assumptions, where available. 



5.3.8 Measuring all of the parameters simultaneously 

Measurements for all of the cosmological parameters are listed in 
Table 17] For each parameter x, we consider both the cases where 
the varied parameters are the vanilla parameters with a minimal ex- 
tension to include the parameter concerned, and the measurement 
in the 'All float' model, when all parameters are simultaneously 
fitted. For the 'All float' model, we compare the result using the 
CMASS P{k) with that using only the CMASS BAO information. 
Comparison of these measurements shows the additional informa- 
tion coming from the power spectrum, even when allowing for nui- 
sance parameters for the galaxy bias modelling using Eq. (|9). 

The ID posterior distribution of the cosmological parame- 



ters and the 2D contours for pairs of these parameters are shown 
in Fig. [191 The constraints using BAO only are generally weaker 
showing that the shape information in P{k) isn't significantly de- 
graded even if a sophisticated galaxy bias modelling with several 
nuisance parameters is used. 



6 CONCLUSIONS 

We have measured the summed neutrino mass using the power 
spectrum of the SDSS-III CMASS sample. We also investigate how 
the prior knowledge of the neutrino mass affects the measurement 
of other cosmological parameters. We do this by studying both the 
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minimal model, where we have a flat, ACDM background with the 
'extra' parameter considered, and also study the bounds on cosmo- 
logical parameters in a more general case, in which all parameters 
are allowed to vary simultaneously. 

For the neutrino mass constraint, we discuss several factors 
which might affect the final result, including the choice of SN data, 
the treatment of the galaxy modelling, and the cutoff scale of the 
power spectrum for the analysis. Our main results can be sum- 
marised as follows, 

(i) How heavy are the neutrinos? 

When we assume a Ami^CDM cosmology, we find that Emj/ < 
0.340 eV (95 percent CL) using a joint analysis of data from 
WMAP7, SNLS3, the CMASS power spectrum, and the BAO 
measurement of WiggleZ, 6dF and SDSS-II. In the most con- 
servative case, where all the cosmological parameters in Table [T] 
are allowed to vary simultaneously, the constraint is weakened to 
Em^ < 0.821 eV (95 percent CL). 

There are several recently published neutrino mass constraints 
listed in Table[8] The more conservative constraints, which was de- 
rived using a smaller /cmax, or assuming a more general cosmology 
and so on, are also listed in the parentheses. The difference among 
these measurements is due to different datasets used and different 
galaxy modelling adopted. The conservative result of this work, 
T^rriu < 0.821 eV is significantly weaker showing how freedom 
in other parameters degrades ^rriu measurements, while the other 
measurements were done in the Ami^CDM cosmology. 

(ii) How many neutrino species ? 

The number of neutrino species is measured to be 4.308 ± 0.794 
in the AA^efrmi^CDM. In this case, the standard value A^eff = 3.046 
is only allowed at the 1.6cr level. In the 'All float' model, the con- 
straint is degraded to A'eff = 4.4121q g^g, but the standard value is 
still only accepted at the 1.57cr level. 

Interestingly, recently several groups have also measured 
A'eff > 3.046, including the Atacama Cosmology Telescope 
(ACT) (Dunkleve t aDl201lh and the South Pole Telescope (SPT) 
(IKeisler e t al. 20JjF teams. This allows the possibility of the exis- 
tence of_one_addkiona^ species, i.e., the sterile neutrino 
rg.g. JCiuffoli et al.ll2012h . 

(iii) Is w = —1? 

If one fits a constant w to the joint dataset, then w — —1.068 db 
0.072 (^CDM) and w = -1.067 zb 0.121 (All float). For the CPL 
parametrisation, which allows w to vary linearly in the scale factor 
a,w — — 1 is still well within the 95 percent CL limit. 

The dynamics of w can be further probed by using a more 
general parameterisatio n, or a non - paran ietric approach to re- 
construct w[a). In fact. Izhao et alj (120 12h recently performed a 
non-parametric reconstruction of w{a) using a similar dataset and 
they found a signal for the dynamics of dark energy at > 2cr level. 
In principle, this is not necessarily inconsistent with our result here 
simply because the CPL parametrisation we used in t his w ork is 
not general enough to cover the behaviour of w{a) that'Zhao et alj 
12012) found. Actually, the CPL parametrisation provides a much 
worse fit to the data (the is only reduced by 0.5 compared to the 
ACDM model) than the lZhao et aP (120121) result, which reduces 
the by 6 with 3 effective number of degrees of freedom. 



(iv) Is the universe close to flat? 

The answer is yes based on this analysis. The conclusion that 
= is robust - the measured Ok is very consistent with zero 
in all the cases we considered. 

(v) Is there a running in the primordial power spectrum ? 

The answer is no. We confirm that as is consistent with zero at 
the 95 percent CL in all the cases we studied. 

(vi) How much tensor perturbation is allowed? 

The tensor- to- scalar ratio is tightened down to r < 0.198 and 
r < 0.440 in the minimal cosmological model and the 'Full fit' 
model respectively. This rules out the 0^ inflation model with the 
e-fold number A^ < 60 to more than 3a. However, the 0^ model 
with A^ > 50 is still allowed. 

All the above issues can be further investigated using the ac- 
cumulating BOSS data in the next few years. 
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